arXiv:1503.08568vl [nucl-ex] 30 Mar 2015 


Spectroscopy of 28 Na: shell evolution toward the drip line 

A. Lepailleur, 1 K. Wimmer, 2,3,4 A. Mutscliler, 5 O. Sorlin, 1 V. Bader, 6 C. Bancroft, 3 D. Barofsky, 3 
B. Bastin, 1 T. Baugher, 6 D. Bazin, 6 V. Bildstein, 7 C. Borcea, 8 R. Borcea, 8 B. A. Brown, 6 L. Caceres, 1 A. 
Gade, 6 L. Gaudefroy, 9 S. Grevy, 10 G. F. Grinyer, 1 H. Iwasaki, 6 E. Khan, 5 T. Kroll, 11 C. Langer, 4 
A. Lemasson, 1,6 O. Llidoo, 1 J. Lloyd, 3 F. Negoita, 8 F. de Oliveira Santos, 1 G. Perdikakis, 3,4 F. 
Recchia, 6 T. Redpath, 3 T. Roger, 1 F. Rotaru, 8 S. Saenz, 3 M.-G. Saint-Laurent, 1 D. Smalley, 4 D. Sohler, 12 
M. Stanoiu, 8 S. R. Stroberg, 6 J.C. Thomas, 1 M. Vandebrouck, 5 D. Weisshaar, 6 and A. Westerberg 3 
1 Grand Accelerateur National d’lons Lourds (GANIL), 

CEA/DSM - CNRS/IN2P3, B. P. 55027, F-14076 Caen Cedex 5, France 
2 Department of Physics, The University of Tokyo, Hongo, Bunkyo-ku, Tokyo 113-0033, Japan 
1 Department of Physics, Central Michigan University, Mt. Pleasant, Michigan 48859, USA 
4 National Superconducting Cyclotron Laboratory, Michigan State University, East Lansing, Michigan 48824, USA 
5 Institut de Physique Nucleaire, IN2P3-CNRS, F-91406 Orsay Cedex, France 
6 Department of Physics and Astronomy and National Superconducting Cyclotron Laboratory, 

Michigan State University, East Lansing, Michigan, 48824-1321, USA 
7 Department of Physics, University of Guelph, Guelph, ON NIG 2W1, Canada 
8 IFIN-HH, P. O. Box MG-6, 76900 Bucharest-Magurele, Romania 
9 CEA, DAM, DIF, F-91297 Arpajon, France 
10 Centre d‘Etudes Nucleaires de Bordeaux Gradignan-UMR 5797, CNRS/IN2P3, 

Universite de Bordeaux 1, Chemin du Solarium, BP 120, 33175 Gradignan, France 
11 Institut fur Kernphysik, Technische Universitdt Darmstadt, 64289 Darmstadt, Germany 
12 Institute of Nuclear Research of the Hungarian Academy for Sciences, P.O. Box 51, Debrecen, H-4001, Hungary 

Excited states in 28 Na have been studied using the /3-decay of implanted 2s Ne ions at GANIL/LISE 
as well as the in-beam 7 -ray spectroscopy at the NSCL/S800 facility. New states of positive 
(J 7 r =3,4 + ) and negative (J 7 r =l-5“) parity are proposed. The former arise from the coupling between 
(M5/2 protons and a 0d3/2 neutron, while the latter are due to couplings with IP3/2 or 0f7/2 neu¬ 
trons. While the relative energies between the J 7 r =l-4 + states are well reproduced with the USDA 
interaction in the N=17 isotones, a progressive shift in the ground state binding energy (by about 
500 keV) is observed between 26 F and 30 Al. This points to a possible change in the proton-neutron 
0 d 5 / 2 - 0 d 3/2 effective interaction when moving from stability to the drip line. The presence of 11^=1- 
4“ negative parity states around 1.5 MeV as well as of a candidate for a state around 2.5 

MeV give further support to the collapse of the N=20 gap and to the inversion between the 0 f 7/2 
and IP 3/2 levels below Z=12. These features are discussed in the framework of Shell Model and 
EDF calculations, leading to predicted negative parity states in the low energy spectra of the 26 F 
and 25 0 nuclei. 

PACS numbers: 21.60.Cs, 23.20.Lv, 27.30.+t 


I. INTRODUCTION 


The first disappearance of a magic shell was proposed 
more than 30 years ago for the neutron magic number 
N=20. This discovery arose from the combined works 
on atomic masses (Tj, nuclear radii [2] and nuclear spec¬ 
tra 13 S3 of nuclei around 32 Mg. This discovery was 
later confirmed by complementary measurements on the 
reduced transition probability values B(E2;0 + -A 2 + ) [3, 
on the quadrupole and magnetic moments [6], as well as 
on the neutron knock-out cross sections [7], to quote a 
few. Theoretically the works of Refs [B3-ETO] described this 
onset of collectivity at N=20 as due to the combination 
of a shell gap reduction and the excitations of particles 
from the normally occupied orbital to the first orbital of 
the upper shell. These excitations lead to a significant 
increase of correlations, eventually bringing the intruder 
configuration from the upper shell below the normal con¬ 
figuration. The nuclei for which the ordering of the in¬ 


truder and normal configurations is inverted belong to 
the so-called ’Island of Inversion’ :T0|. Their configura¬ 
tions are often strongly mixed. 

The shell evolution of neutron orbits as a function 
of the proton number is illustrated in Fig. [l] for the 
N=20 isotones. This figure has been elaborated with the 
monopole proton-neutron interactions obtained in the 
shell model approach of Ref. m that include the few 
experimental observations made at that time. At Z= 8 , 
the neutron 0 d 3 / 2 , IP 3/2 and 0 f 7 /2 orbits are unbound, 
the N=20 is weaker than the N=16 gap, and the order¬ 
ing of the IP 3/2 and 0 f 7 / 2 orbits is reversed compared to 
what is observed in the valley of stability. While filling 
the proton 0d 5 / 2 orbit from Z =8 ( 28 0) to Z=14 ( 34 Si), 
the N=20 gap increases as the N=16 gap is reduced. Af¬ 
ter Z=14, the filling of the ls^ and 0d 3 / 2 proton or¬ 
bits keeps the N=20 gap unchanged between Z=14 and 
Z=20. The formation of the N=20 shell gap, as well as 
the inversion between the neutron 0 f 7 /2 and IP 3/2 or¬ 
bitals are profound structural changes in the shells that 
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FIG. 1: (color on line) Effective neutron single-particle en¬ 
ergies (ESPE) of the N=20 isotones between Z =8 and Z=20 
based on the work of Ref. nu. The slope of the lines corre¬ 
sponds to the strength of the proton-neutron monopole in¬ 
teractions. The red line, which has the largest slope, shows 
the effect of the strong 65 / 2 - 63/2 proton-neutron monopole 
interaction when the proton 65/2 is filled. 
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are likely caused by the hierarchy between the proton- 
neutron monopole interactions involved: the 0d 5 / 2 -0d 3 / 2 
interaction is much larger than the others involved such 
as the 0ds/ 2 -0f 7 / 2 and 0 d 5 / 2 -lp 3 / 2 ones, as described in 
Refs. [TTl - iToj . When the lp 3 / 2 orbit becomes unbound a 
further reduction of the monopole interactions involving 
this \crw-l orbital may be present. 

The 0d 5 / 2 -0d 3 / 2 , 0d 5 / 2 -lp 3 / 2 and 0d 5 / 2 -0f 7 / 2 monopole 
interactions are so far poorly constrained by experimen¬ 
tal data. To achieve this goal, the evolution of the 0d 3 / 2 , 
lp 3 / 2 and 0f 7 / 2 neutron single-particle energies should 
be determined from the Si to O isotones. However the 
N=20 nuclei lying in the island of inversion are deformed. 
Consequently their configurations are strongly mixed and 
their spherical single-particle energies can rarely be de¬ 
termined. A further difficulty arises from the fact that 
the N=20 nuclei with Z=10 and Z=8 are weakly bound 
or unbound. They are difficult to produce in experiments 
and their study would often require the treatment of the 
continuum. In order to study these interactions and their 
behavior toward the neutron drip line, we take advan¬ 
tage of the fact that large N=14 and N=16 shell gaps are 
present around Z=8-10 [HH20] . where the N=20 gap is 
weakened. The N=17 nuclei can therefore be described 
mainly with a single neutron in the d 3 / 2 , p 3 / 2 or f 7 / 2 
orbit, without a large mixing with neighboring shells. 

Recent experimental observations pTNBTj have brought 
remarkable credit to the theoretical description of Fig [I] 
around (N=16, Z=8) suggesting an inversion between the 
lp 3 / 2 and 0f 7 / 2 orbits in the N=15 and N=17 Ne iso¬ 
topes. These experiments have also demonstrated that 
the N=20 gap is small. In addition, the works of ITS1-I2()j 


have proven that the N=16 shell gap is large at Z=8, 
while Ref. [25] suggests that the N=16 gap somehow 
persists at Z=10 in 26 Ne, as witnessed by its vibrational 
behavior. 

The best way to extract information about the afore¬ 
mentioned monopole and multipole parts of the nuclear 
interaction is therefore to study odd-odd nuclei in the 
N=17 isotones. The coupling of protons and neutrons in 
the d 5 / 2 and d 3 / 2 orbits will lead to J 77 =l-4+ states in 
the f 3 Ali 7 , ?fNai 7 and g 6 Fi 7 isotones that span from near 
stability to the neutron drip line. In a similar manner, 
the coupling of the protons in the d 5 / 2 orbit with neu¬ 
trons in the lp 3 / 2 (0f 7 / 2 ) orbit leads to J w =l-4 _ (J”’=l- 
6 _ ) negative parity states in the same nuclei. However 
such studies are tedious as several experimental meth¬ 
ods are often required to produce all the states of these 
multiplets for nuclei that are not so easily produced at 
radioactive ion beam facilities. The first J=l-4 positive 
parity states in 30 A1 have been obtained recently using 
the Gammasphere array [26] . Negative parity states are 
proposed from 2.29 MeV on, but they were assumed to 
originate from the neutron f 7 / 2 orbital. In 26 F three dif¬ 
ferent experimental techniques were required to study the 
bound J 77 =2,4 + states [25, .27] and the unbound J 7r =3 + 
state [25] - So far there is no evidence of negative par¬ 
ity states in 26 F below the neutron emission threshold 
of 1.070(62) MeV [25] ■ As for the 28 Na nucleus, there 
are only candidates for the J 77 = l) 1 " and J 77 = 2f states 
that were derived from a previous /3-decay study of 28 Ne 
[29] . The present study aims at determining the energy 
of the missing J 77 = and J 77 = 4f states and providing 
information on the presence of negative parity states in 
28 Na to confirm the lowering of the lp 3 / 2 and 0f 7 / 2 orbits 
toward Z=8. To achieve these goals two complementary 
experimental techniques were required. We first repeated 
the /3-decay experiment of 28 Ne at the GANIL/LISE fa¬ 
cility with a larger statistics as compared to 23 and we 
secondly used the in-beam spectroscopy technique at the 
NSCL/MSU facility to detect the q-rays of 28 Na in the 
GRETINA Ge detector array [3U] produced in the neu¬ 
tron and proton removal reactions from 31,32 Mg. 

II. /3-DECAY EXPERIMENT 
A. Experimental technique 

The 28 Ne nuclei were produced at the GANIL facility 
through the fragmentation of a 77.6 A • MeV 36 S 16+ beam 
with an intensity of 2 /iAe on a 237 mg/cm 2 Be target. 
Nuclei of interest were selected by the LISE m spec¬ 
trometer, in which a wedge-shaped Be degrader of 1066 
/im was inserted at the intermediate focal plane. They 
were identified from their energy losses in a stack of three 
500 ;tm Si detectors and from their time-of-flight refer¬ 
enced to the cyclotron radio frequency. The spectrometer 
set-up was optimized for 26 F, but a fraction of the 28 Ne 
ions were also transmitted at a rate of 10 pps, corre- 
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FIG. 2: Identification of the nuclei produced in the experi¬ 
ment through their energy loss and time of flight. 


f onding to 22% of the implanted nuclei shown in Fig. 

A total of 3.24 xlO 6 28 Ne nuclei were implanted in a 
1 mm-thick double-sided Si stripped detector (DSSSD) 
composed of 256 pixels (16 strips in both the X and Y 
directions) of 3x3 mm 2 located at the final focal point of 
LISE, ,5-particles were detected in the same strip of the 
DSSSD as the precursor nucleus 28 Ne. 

Four Ge detectors of the EXOGAM array 132] sur¬ 
rounded the DSSSD to detect the 7 -rays with an effi¬ 
ciency of 6.5(3)% at 1 MeV. A total /3-efficiency of 60(1)% 
has been determined from the comparison of the inten¬ 
sity of a given 7 -ray belonging to the decay of 28 Ne gated 
or not on a /3-ray. A Si-Li detector of 5 mm thickness 
was placed downstream of the DSSSD for three purposes: 
i) to ensure that the 28 Ne ions were implanted in the 
DSSSD and not passing through it, ii) to discriminate 
very light particles in the beam that pass through the 
telescopes and iii) to determine the /3 energy threshold 
in each strip using the coincidence with the ^-particles 
detected in Si-Li detector. 


B. Beta-decay scheme of 28 Ne 

The upper (bottom) spectrum of Fig. [3] displays the 
7 -ray energy spectrum obtained with a /3-particle corre¬ 
lated in space and time up to 40 ms (between 120 and 
1000 ms) after the implantation of a 28 Ne precursor. Ow¬ 
ing to the short lifetime of 28 Ne (see below), the transi¬ 
tions from Fig. [3] belonging to its /?-decay are more in¬ 
tense in the upper (blue) spectrum, while those caused 
by daughter decays or other implanted nuclei are domi¬ 
nating the bottom (red) spectrum. Consistent half-lives 
values ofT 1/2 =18.2(5) ms and 18.6 (2) ms are found from 
the /3-decay time spectrum of 28 Ne gated on the 863 keV 
and 2063 transitions seen in Fig. [3j respectively. These 
values are also in accordance with T 1 / 2 ( 28 Ne) = 18.4 (5) 
ms derived in Ref. [20] . 

The level scheme shown in the left part of Fig. [ 8 ] 
is established from /3-gated 7-7 coincidences following 
the implantation of a 28 Ne nucleus. Based on their di¬ 
rect /3-decay feeding from the 0 + ground state of 28 Ne 
and Gamow-Teller /3-decay selection rules, we propose a 


TABLE I: List of states, 7 -ray energies and relative intensities 
(normalized to 100 28 Ne decays) observed in the /3-decay of 
28 Ne to 28 Na. 


Ei (keV) 

E{ (keV) 

E . 7 (keV) 

J rel (%) 

55 a 

0 

55 

- 

691(1) 

55 

636(1) 

0.3(1) 

1131(1) 

0 

1131(1) 

0 . 2 ( 1 ) 


55 

1076(1) 

1 . 0 ( 1 ) 

1255(1) 

0 

1255(1) 

2.7(2) 


55 

1200 ( 1 ) 

0.5(1) 


691 

564(1) 

0.3(1) 

1932(1) 

0 

1932(1) 

0.5(1) 


55 

1877(2) 

0.4(1) 

2118(1) 

0 

2118(2) 

0 . 8 ( 1 ) 


55 

2063(1) 

14.2(12) 


55 

863(1) 

3.6(3) 

2714(1) 

0 

2714(2) 

0.9(1) 


55 

2659(1) 

0.9(1) 


1131 

1583(1) 

0.7(1) 


1932 

782(1) 

1 . 1 ( 1 ) 


2118 

596(1) 

0 . 6 ( 1 ) 


“taken from Ref. m 


r= 1 + spin parity value to the states at 2714, 2118 keV 
and to the ground state of 28 Na. The /3 feedings to the 
28 Na ground state as well as to unbound states (leading 
to the /3-delayed neutron emission) are derived from the 
intensities of the observed lines seen in Fig. [^populating 
excited states in the 28 Mg and 2 'Mg nuclei, respectively. 
They agree with the experimental feedings determined in 
Ref. [29] for most of the states, including for unbound 
states. As shown in Fig[ 8 j these /3 feedings compare rea¬ 
sonably well with shell model calculations. However, as 
for the /3-decay of 26 Ne [323 the feeding of only three of 
the four predicted 1 + states is observed in the /3-decay 
of 28 Ne below the neutron separation energy. 

With the exception of the 3231 keV and 3457 keV tran¬ 
sitions, all the transitions belonging to the 28 Ne decay to 
28 Na observed in j22] are observed here. With /3 feed¬ 
ings of ~ 2 % for the 3231 keV and 3457 keV transitions 
given in Ref. [29], they should have been seen in our data 
with a confidence level of 10 er. As the spatial correlation 
was not used in Ref. [29], it is possible that these 7 rays 
were wrongly assigned to the 28 Ne decay. We observe 
two new 7 rays at 564(1) keV and 636(1) keV of similar 
weak intensity that are in mutual coincidence (see Fig. 
[4]) as well as in coincidence with the 863(1) keV 7 rays 
de-exciting the 2118(1) keV level in 28 Na. The summed 
energy of these two 7 transitions, 1200(2) keV, matches 
the energy of the 1200(1) keV transition coming from the 
decay of the 1255(1) keV to 55.2(5) keV states in 28 Na. 
We therefore propose a new level in 28 Na, the placement 
of which (691(2) keV) was derived using the information 
obtained in a second experiment described below. The 
list of states, 7 -ray energies and intensities observed in 
the /3-decay experiment is given in Table [I] 

To propose spin parity assignments to the identified 
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FIG. 3: (color on line): /3-gated 7 -ray spectra obtained up to 40 ms (blue, upper) and between 120 ms and 1000 ms (red, 
bottom) after the implantation of a 28 Ne nucleus. The red histogram has been scaled down to see the two spectra on top of 
each others. Identified 7 -rays are shown with different symbols. The new ones at 564 (1) and 636 (1) keV attributed to the 
decay of 2S Ne are visible on the left hand side of the spectrum. 



FIG. 4: (color on line) Left: /3-gated 7 -ray spectrum obtained 
up to 60 ms (blue) and between 120 ms and 1000 ms (red) 
after the implantation of a 28 Ne nucleus. Right: /3-gated 7 -ray 
spectrum in coincidence with the 564 (1) keV transition. 


the J 7 r =2 + at 55 keV with a 7 -ray of 636 keV rather than 
through an E2 transition to the J 7 r =l + ground state. We 
propose spin parity assignments ^=?>2 and J 7 r = 2 j to 
the states at 1131 keV and 1255 keV from the compari¬ 
son to shell model calculations as well as from their feed¬ 
ing and decay branching ratios in E2 and Ml transitions. 
The 1932 keV state likely has a J’ r =2 + configuration as 
it decays equally to the 1 + ground state and to the 2 + 
state at 55 keV. Most of the states predicted by the shell 
model calculations using the USDA interaction have their 
experimental equivalent up to 1255 keV. 


III. IN-BEAM 7 - RAY SPECTROSCOPY 


states, we start with the fact the ground state of 28 Na 
has J 7 r =l + , which is deduced from its large direct feeding 
from the 0 + ground state of 28 Ne. The excited state at 55 
keV is likely to have J^O,2 + as a larger spin difference 
between this 55 keV state and the ground state would 
have resulted in a long-lived isomer transition neither 
observed in the /?-decay part of our work nor in Ref. [25] . 
Between the J 7 r =0,2 + candidates, the J 7 r =2 + value seems 
the most reasonable, considering that the USDA [32] pre¬ 
diction gives the first 2 + state around 100 keV and the 
first 0+ state at a much higher energy of about 2 MeV. 
Based on the comparison to shell model calculations and 
their respective feeding from higher energy states and 
decay branches to the 55 keV or ground state, the state 
at 691 keV is a good candidate for J 7 r =3 + . Indeed a 
J 7r =3+ state mainly decays through an Ml transition to 


A. Experimental technique 

In-beam 7 -ray spectroscopy of the neutron-rich 28 Na 
isotope was performed in a second experiment at the 
Coupled Cyclotron Facility of NSCL at Michigan State 
University. The 28 Na nuclei were produced in the sec¬ 
ondary fragmentation reactions from 31 Mg and 32 Mg 
beams impinging at about 95 MeV/u on a 375 mg/cm 2 
°Be target. In total about 7.7 ■ 10 5 7 — 28 Na coincidences 
were recorded in the two settings together. The outgoing 
particles were identified based on the time-of-flight and 
energy-loss measurements using the focal-plane detection 
system of the S800 spectrograph [35]. Trajectories of re¬ 
coiling ions were tracked in the S800 using two sets of 
Cathode Readout Drift Chambers measuring their po¬ 
sition and angle values at the focal plane allowing the 
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reconstruction of their velocity (/3 ~ 0.41) and position 
at the target. The target was surrounded by seven mod¬ 
ules of the GRETINA array [SHI, each module consisting 
in four HPGe segmented crystals. Four modules were 
placed at the most forward angles around 58° and three 
around 90° to detect 7 rays induced by the de-excitation 
of the nuclei in-flight with an efficiency of about 6 % at 1 
MeV. The 7 -ray detection threshold has been lowered to 
about 50 keV in order to be able to observe the low-lying 
transition from the first excited state in 28 Na at 55 keV. 
The threshold and energy resolutions were modeled using 
the GEANT4 [3S] simulation of the GRETINA setup [37]. 
The 28 Na velocity vector and position were used to apply 
an event-by-event Doppler correction to the 7 -ray spec¬ 
trum of the 28 Na nuclei shown in Fig. [ 5 J A 7 -ray energy 
resolution (FWHM) of about 2% at 1 MeV was achieved. 
The list of 7 -ray energies and intensities observed in this 
experiment is given in Table |H| 

TABLE II: List of states, 7 -ray energies and intensities ob¬ 
served in the in-beam experiment. Systematic uncertainties 
a(E) of about 3 keV should be considered on the proposed 
7 -ray energies. 


E-, (keV) 

E t (keV) 

E 1 (keV) 

I lel (%)“ 

55 

0 

55 

380(50) 

688 

55 

633 

100 ( 2 ) 

968 

55 

912 

7(1) 


688 

277 

52(3) 

1131 

0 

1131 

1.9(6) 


55 

1075 

16.6(8) 

1233 

0 

1233 

1.9(6) 


55 

1177 

29(1) 


688 

542 

8.3(6) 

1255 

0 

1255 

15.7(8) 

1353 

55 

1298 

33(1) 

1481 

0 

1481 

8 . 2 ( 8 ) b 

1636 

688 

948 

11 . 6 ( 8 ) 


1233 

403 

16.7(7) 


1353 

282 

10 ( 2 ) 

1749 

55 

1694 

5.2(7) 

1792 

0 

1792 

6.9(9) 

1929 

0 

1929 

8.5(9) c 

2121 

55 

2066 

8 . 2 ( 8 ) c 

2378 

0 

2378 

5.2(8) d 

2493 

965 

1527 

4(1) 


1636 

858 

14.8(8) 

2605 


2605 

10 ( 1 ) d 

2650 


2650 

7(1) d 

2874 


2874 

4(1) d 



489 

5.7(6) d 


“normalized to the 633 keV transition 
^placement uncertain 
“placement based on ft decay 
d unplaced 

The systematic uncertainty for the intensity of the 55 
keV transition is particularly large because of the vicin¬ 
ity of the threshold. Additionally, this state exhibits a 
relatively long lifetime as deduced from its low energy 



FIG. 5: (color on line) Doppler-reconstructed 7 -ray energy 
spectrum obtained from in-beam 7 -decay spectroscopy of 
28 Na. The solid green line is the response obtained from the 
GEANT4 simulation if individual transitions (blue lines) and 
a continuous background (orange line) are considered. The 
inset shows the high energy part of the spectrum with in par¬ 
ticular the doublet of transitions at 2605 and 2650 keV. 


tail. It follows that the Doppler shifted energy measured 
in GRETINA lies partially below the detection thresh¬ 
old. The lifetime of the 2 + state has been simulated, 
and the resulting response function has been fitted to 
the experimental spectrum. This way a lower limit for 
the lifetime of t > 1 ns is obtained. An additional con¬ 
straint on the lifetime can be obtained from the fact that 
100 % of the 633 keV 7 -ray yield has to proceed through 
the 55 keV transition. Consistent yields in the 633-55 
keV coincidence are obtained for simulated lifetimes of 
r = 1.4(4) ns. Shell model calculations using the USDA 
effective interaction predict the 2 + state at 182 keV. 
Using the shell model reduced transition probabilities, 
B(E2) and B(M1) values, for the decay to the ground 
state, and the experimental transition energy of 55 keV, 
a theoretical lifetime of Ttheo = 1-42 ns is obtained. This 
validates our assumption of a relatively long lifetime of 
the 2 + state at 55 keV. 

Due to the high detection efficiency of the GRETINA 
array it was possible to construct 7 — 7 coincidence spec¬ 
tra by gating on several transitions. Fig. [ 6 ] shows the 
spectrum observed in coincidence with the 55 keV tran¬ 
sition. This spectrum shows which of the observed 7 
transitions are directly or indirectly populating the first 
excited state at 55 keV. 


B. Level scheme of *Na 

The level scheme of 28 Na as well as tentative spin as¬ 
signments derived from this part of the experiment are 
based on the single 7 -ray energy spectrum of Fig. § 7-7 
coincidence spectra of Fig. [ 6 ] and [7] as well as 7 branching 
ratios from Table [IT] Two almost overlapping 7 transi¬ 
tions are present at 277 and 282 keV. The 633(3) keV 
7 -ray is compatible with the 636(1) keV 7 -ray observed 
in the /3-decay experiment. As the 564(1) keV 7 -ray 
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FIG. 6 : (color on line) Doppler-reconstructed 7 - 1 'ay energy 
spectrum gated on the 55 keV transition from the first excited 
state in 28 Na in two different energy ranges. The blue line 
(top spectrum in each pannel) shows the total spectrum for 
comparison. While most of the transitions are observed in 
coincidence with the 55 keV 7 ray, the transitions of 1255 
and 1481 keV are for instance clearly missing in the 7 -gated 
spectrum. 


found in the /3-decay experiment is not observed here, 
it is placed above the 636(1) keV transition. The 277(1) 
keV 7 -transition is in coincidence with the 636 keV one. 
As the intensity of the latter is larger, it is placed below 
the 277 keV transition, establishing a new state at 968(3) 
keV as seen in Fig. [ 8 ] A tentative 3 7 I ’=4 + assignment is 
proposed for this level as it mainly decays by a 277 keV 7 - 
ray to the previously assigned 3 + at 691 keV, very weakly 
to the 55 keV 2 + state and not to the ground state with 
J 7 r =l + . With this newly proposed J 7 r =4 + state, all pos¬ 
itive parity states but the J 7 r =0 + that are predicted by 
the shell model calculations below 2.2 MeV are observed 
experimentally. As candidates for 3 7 r =0 + state are pro¬ 
posed in the following paragraph most of the other states 
populated in this reaction and shown in the right part of 
Fig-® are proposed to be candidates for intruder states, 
arising from the neutron fp shells. 

A tentative 3 7 r =2 _ assignment is proposed to the 1233 
keV level from the fact that it decays to the J 7 r =l + 
ground state and to the 3 7 r =2+ state at 55 keV. A J’ r =3 _ 
assignment is proposed to the 1353 keV level on the ba¬ 
sis of its sole decay to the J 7 r =2 + state through a 1298 
keV transition. Owing to the fact that the 1481 KeV 
state exclusively decays to the ground state, its spin as- 
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FIG. 7: (color on line) 7 — 7 coincidence spectra of 28 Na 
obtained during the in-beam experiment when gating on sev¬ 
eral 7 -rays (from top to bottom 633, 403, 858 and (633,1298) 
keV). 


signment could be 3=0-2, providing a good candidate to 
the missing J’ r =0 + state. A new level is firmly estab¬ 
lished at 1636 (2) keV from its observed three 7 decay 
branches 282 + 1298 + 55, 403 + 1177 and 948 + 636 
+55 keV. As this level decays to the 3 7 r =3“, 3’ r =2~ and 
3 7 r = 3 + states and not to the 3=1,2 positive parity states 
at lower energy, it is likely to have 3^=4“. Two lev¬ 
els with tentative spin 3 7 r =l~ and 3 Tr =2 _ are proposed 
at 1749 and 1792 keV, respectively. The level at 1792 
keV could as well be a candidate for the 3 7 r =0 + state 
predicted at a similar energy by the shell model calcu¬ 
lations. A 3 7 r =5 _ level is proposed at 2493 keV from 
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its 858 and 1527 keV 7 branches to the previously as¬ 
signed J n =4~ and J”'=4 + states at 1636 and 968 keV, 
respectively. Though less probable, a J 7 I '=3 _ assignment 
cannot be ruled out. Other high energy 7 -rays are ob¬ 
served at 2378, 2605, 2650 and 2874 keV. However, their 
placement in the level scheme is uncertain by 55 keV as 
there is not enough statistics to ensure that a coincidence 
with the 55 keV 7 -ray is present or not. Finally, the 489 
keV 7 -ray could not be placed without ambiguity in the 
level scheme as well. 


IV. DISCUSSION 

A. Positive parity states along the N=17 isotones 

The level scheme of 28 Na shown in Fig. | 8 ]has been ob¬ 
tained by combining the results from the /3-decay (left 
part) and the in-beam 7 spectroscopy (right part) exper¬ 
iments. The middle part of the spectrum displays tran¬ 
sitions that were common to the two experiments.The 
two new levels with configurations J 7r = 3]~ and 4]" com¬ 
plete the quadruplet of J 1 I '=l-4 + states resulting from the 
7 rd 5 / 2 8 ^d 3 / 2 coupling in 28 Na. Added to the recently 
discovered J ir =l-4 + states at low energy in the N=17 
isotones of 26 F [23 [27, .23] and 30 A1 [26] , a systematics 
of their binding energies as a function of the proton-to- 
neutron binding energy asymmetry can be obtained and 
is compared to shell model calculations in Fig. [9] The 
energy difference between proton and neutron separation 
energies S p -S„ ranges from 15 MeV in the close-to-drip- 
line nucleus 26 F to 7 MeV in 30 Al. For each nucleus 
the calculated binding energy of the ground state (using 
the USDA interaction) is taken as the reference value for 
Fig. [9] 

The J 7 r =l-4 + states in the 26 F nucleus can be de¬ 
scribed as a proton 7 rd, 5/2 coupled to a neutron in the 
1 ^ 3/2 on top of a 24 0 core nucleus. In this particle- 
particle coupling scheme, the multiplet of states in 26 F 
displays an upward pointing parabola in binding energy 
value as a function of J as seen in Fig. [9] In the present 
case the amplitude of the parabola scales to a first order 
with the strength of the residual interaction that splits 
the different components of the J 7 r =l-4 + multiplet. With 
the exception of the J 7 r =3 + state that is unbound and 
therefore may need a specific treatment, the ordering of 
the states in 26 F is well reproduced by the USDA inter¬ 
action but their calculated binding energies are too large 
by about 200 keV. Even if the J 7 I '=l-4 + states in 30 Al 
do not have a pure configuration, they can be viewed as 
hole-particle coupling 8 1 M 3/2 states with re¬ 

spect to the full occupancy of the proton ds /2 orbit. This 
leads to the downward pointing parabola as a function 
of J shown in Fig. [9] The amplitude and shape of the 
parabola is extremely well reproduced using the USDA 
interaction. However, the calculated absolute binding en¬ 
ergies are this time smaller than the experimental values 
by about 300 keV. The pattern of the binding energy of 


the J ir =l-4 + states in 28 Na is intermediate between 26 F 
and 30 Al owing to the fact that the ds /2 orbit is only 
half-filled. The spectrum of 28 Na is well reproduced by 
the calculations. It is worth to notice here that a similar 
shift in binding energy is observed along the N=17 iso¬ 
tones when using the USDB interaction. However such 
calculation proved to be less precise, predicting the bind¬ 
ing energy of the J 7 r =4+ state wrong by about 600 keV 
in 26 F and 400 keV in 30 Al and a J 7 r =2 + ground state of 
28 Na instead of J=l + . 

Gathering all results along the N=17 isotonic chain, a 
systematic deviation between experimental and theoret¬ 
ical binding energies is observed: as compared to experi¬ 
mental values the 26 F is over bound, the 28 Na is perfectly 
well reproduced, while the 30 Al is under bound. A to¬ 
tal shift in energy of about 500 keV is found between 
these three isotones. Being a global shift in energy of the 
whole multiplet, this change in binding energy is likely 
carried by the monopole part of the 7 rd 5/2 <8 ZM 3/2 nu¬ 
clear interaction rather than by multipoles. We propose 
three reasons that may explain this systematic deviation. 
First, on the experimental side, the atomic mass of 26 F 
derived in [38], that has been corrected in Ref. [25] to 
account for the possible contamination from the 4 + iso¬ 
meric state, may not be correct. This would account for 
the shift observed in 26 F but not in 30 Al. Second, the 
shell model interaction of [34] does not use an explicit 
isospin dependence of the Coulomb energy contribution 
to the binding energy of the nuclei as it should probably 
be done to account for the change of nuclear radii for a 
wide range of isotopes. This part is therefore taken only 
implicitly in the fitting procedure to determine TBME 
from the experimental data, a feature than can cause 
some bias in the calculation of the binding energies if the 
suitable data to constrain this isospin dependent part are 
not available experimentally or if they are not included 
in the fitting procedure. Third, the scaling law of the 
monopole interaction in A ” 1 / 3 may not be appropriate, 
and more realistic interactions should be developed for 
nuclei spanning over a large proton-to-neutron binding 
energy asymmetry. 


B. Negative parity states in the N=17 isotones 

In the odd-even nuclei N=17 isotones, negative par¬ 
ity states 7/2“ and 3/2“ originating from the O/7/2 
and lp 3 /2 orbits have been populated at about 3 MeV 
through (d,p) reactions: in 33 S, the first 7/2“ and 3/2“ 
states lie at 2.934 and 3.220 MeV with C 2 S values of 
0.53 and 0.87, respectively [31]. In 31 Si, C 2 S values for 
the first 7/2 - and 3/2~ at 3.133 and 3.535 MeV are 
0.6 and 0.4, respectively [39] ■ Recently, negative par¬ 
ity states have been observed in 29 Mg and 27 Ne. Their 
excitation energies are much lower, and their ordering 
are reversed compared to nuclei in the valley of stability. 
States with L =1 and L=3 have been populated in 29 Mg 
at energies of 1096 and 1432 keV using the one neutron 
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FIG. 8 : (color on line) The very right part displays shell model predictions using the USDA (WBA-M) interaction for the 
positive (negative) parity states. The rest shows the experimental level scheme of 28 Na derived by adding information from 
the /3-decay of 28 Ne (left) and from the in-beam gamma-ray spectroscopy (right). The middle part displays 7 -lines that are 
observed in the two experiments. The energies and uncertainties of the levels are derived from the /3-decay when possible (Table 
[I|, from the in-beam experiment (Table [TlJ) , or from a combination of the two experiments otherwise. 



FIG. 9: (color on line) Comparison between the experi¬ 
mental binding energies of the J n = 1 + - 4 + states in 
30 A1 [26[ |40] (red), 28 Na (29[ [40] and present work (green), 
and in 26 F |25l 12711281I38| (blue) and the USDA shell model 
predictions (in black, with the g.s. value written for each nu¬ 
cleus). The reference value for each nucleus is the calculated 
binding energy of the ground state. Experimental binding en¬ 
ergy uncertainties for the 30 Al, 28 Na and 26 F ground states 
are 14, 10 and about 130 keV, respectively. 


knock-out reaction from 30 Mg [41J . Similar L assign¬ 
ments have been proposed to the states at 0.765 (3/2“, 


L=l) and 1.74 MeV (7/2“, L=3) that are populated in 
the 26 Ne(d,p) 27 Ne reaction with C 2 S values of 0.64(33) 
and 0.35(10), respectively [S2J. Globally as protons are 
removed from the 0d^/2 orbit from 31 Si, the excitation en¬ 
ergy of the two negative parity states 7/2“ and 3/2“ is 
decreasing relatively to the 3/2+ ground state, and their 
ordering is reversed around Z=12. It is derived that the 
N=20 gap between the 0d 3 / 2 and 0/ 7 / 2 is collapsing and 
the traditional N=28 gap between the 0/ 7 / 2 and 1 p 3 / 2 is 
upside down. These features have been qualitatively at¬ 
tributed to the hierarchy of the proton nuclear forces in 
Ref. US]. In the odd-odd nuclei such as 28 Na, multiplet 
of negative parity states J 7r =l _ -4^ and J 7r =2~-5“ are 
expected to be formed by the proton-neutron coupling 
0c? 5 / 2 -lp 3 / 2 and 0d 5 / 2 -0/ 7 / 2 , respectively. The present 
observation of tentatively assigned negative parity states 
J=l,2,3 and 4 in 28 Na at 1749, 1233, 1353 and 1636 keV, 
respectively, confirms the presence of the neutron lp 3 / 2 
orbit at relatively low energy at Z=10 [22] . Actually 
the J=l-4 negative parity states are systematically more 
bound than predictions by about 300 keV, hinting at a 
possibly stronger proton-neutron Od 5 / 2 -lp 3 / 2 interaction 
than calculated. The tentatively assigned 3 lr =5~ state 
at 2496 keV (calculated at 2867 keV) is in accordance 
with the presence of the 0/ 7 / 2 orbit about 1 MeV above 
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the IP 3/2 orbit, as proposed in 27 Ne [22] . 


C. Evolution of nuclear structure toward drip-line 

In order to understand the evolution of the IP 3/2 or¬ 
bit and of the N=20 gap toward the neutron drip line, 
configuration-interaction (Cl) calculations have been car¬ 
ried out with a Hamiltonian called WBA. It is the same 
as the WBP Hamiltonian from [321, but the older USD 
sd-shell Hamiltonian [43] part has been replaced by the 
more recent USDA Hamiltonian [33] that was used in the 
previous section to calculate the energy of positive-parity 
states. The required basis to model the N=17 isotones 
is the full (IsOd) for positive parity states with a core of 
16 O. Moreover it allows for one neutron to be excited to 
the (lpOf) valence space for negative parity states. While 
this basis can be used for low-lying negative-parity states 
in 25 0 , 26 F or 27 Ne, the dimension of the calculation is 
at the limit for 28 Na and too large for 29 Mg and 31 Si. 
Such calculations might be possible in the future. 

The WBA interaction is used to calculate the energies 
of negative-parity states in 27 Ne. As they are too high 
by 0.4 MeV, the single-particle energies of the lp-Of shell 
orbitals have been modified by 0.4 MeV, leading to the 
WBA-M interaction. In Table |III| excitation energies of 
low-lying negative parity states are given for the N=17 
isotones that the model can handle. As used to con¬ 
strain the WBA interaction, the spectroscopy of 27 Ne is 
well reproduced. Negative parity states are present at 
low energy in 28 Na, as found experimentally. The wave 
functions of the first J w =l-4~ are composed by 60-75% 
of a neutron in the lp 3 / 2 and by 20-30% of a neutron 
in the O/ 7/2 orbital, the remaining weak fraction being 
distribution in the lpi / 2 and 0 f 5 / 2 orbitals. All calcu¬ 
lated negative parity states in 26 F lie above the neutron 
emission threshold of S„ = 1.070(62) MeV [25], a fea¬ 
ture that agrees with the observation of only two bound 
excited states of positive parity in this nucleus. In 25 0, 
the 3/2“ is expected to lie only 493 keV above the 3/2+ 
resonance [19]. The three negative-parity states 3/2“, 
7/2“ and 1/2“ in 25 0 have large 24 0+n spectroscopic 
factors. With the proximity in energy between the 3/2 + 
and 3/2“ states in 25 O, 26 O likely contains a significant 
amount of negative parity contribution in its ground and 
first excited states. The reliability of this extrapolation 
far from stability depends on the confirmation of spin as¬ 
signments of negative parity states in 28 Na, on the pos¬ 
sibility to model higher Z isotones in this large valence 
space, as well as on the possibility to observe of negative 
parity states in 26 F in the future. 

The present shell evolution towards the neutron drip 
line is made in an Harmonic Oscillator basis in which 
bound and unbound states are treated on the same foot¬ 
ing. Therefore we look at trends of ESPE in the frame¬ 
work of energy density functional (EDF) calculations 
that do not use a Harmonic Oscillator basis. The single¬ 
particle energies obtained with the Skx [33] and Skxtb 


TABLE III: Excitation energies of some low-lying negative- 
parity states obtained with the WBA-M Hamiltonian are 
compared to experimental values. Note that the spin assign¬ 
ments of the highest energy states in 28 Na are tentative, see 
text and Fig. [8] for details. 


Nucleus 

r 

E* (MeV) 
theory 

E* (MeV) 
experiment 

25 O 

3 / 2 - 

0.493 



1 / 2 - 

1.898 



7/2" 

2.611 


2Ep 

4~ 

1.339 



2~ 

1.384 



1~ 

1.952 



3~ 

2.485 


27 Ne 

3/2" 

0.825 

0.765 [23] [24] 


7/2" 

1.710 

1.74 [22] 


1/2- 

1.834 


2S Na 

2" 

1.552 

1.233 


3- 

1.715 

1.353 


4- 

1.888 

1.636 


1“ 

2.100 

(1.749) 


2" 

2.341 

(1.792) 


3- 

2.821 



O' 

2.836 



5- 

2.867 

2.493 


[35] Skyrme functionals are shown in Fig. 10 for an N=16 
core as a function of proton number. This figure is sim¬ 
ilar to Fig. [TJ the major changes consisting here in the 
down bending of the IP 3/2 orbital when reaching the drip 
line and a weaker reduction of the N=20 gap (between 
0d 3 / 2 and Ofy / 2 orbits). This non linearity in the ESPE 
of the IP 3/2 orbital is due to its halo-like nature that is 
reducing the effective monopole interactions there. 


The ESPE for an unbound state was estimated as fol¬ 
lows. First the central EDF potential is increased by a 
factor that binds the state by 0.2 MeV. Then the ex¬ 
pectation value of the full kinetic and potential energy of 
this state is calculated resulting in a positive value for the 
ESPE that is shown. (More correctly, the properties of 
unbound states should be related to position and shape of 
neutron scattering resonances.) The weak-binding effect 
in the single-particle energies for the low-£ orbitals for 
light nuclei has been discussed by Hamamoto [32] (see 
Figs. 3 and 4 in [12]). It is the same mechanism that 
reduces the lsi/ 2 -0d 5 / 2 gap for the region of Z=2-8 as 
shown by the EDF calculations in Fig. 1 of m, and 
more recently by calculations with a Woods-Saxon po¬ 
tential [35]. Basically the binding energy of the low £ 
orbits is more significantly affected by the proximity of 
continuum as they encounter a much smaller centrifugal 
barrier than higher l orbits. 

Fig. [To] also shows the results of two different calcu¬ 
lations. Skx does not contain a tensor interaction, while 
for Skxtb a tensor interaction has been added that de¬ 
scribes the trends of single-particle energies in heavier 
nuclei 05]. The ESPE for Skx and Skxtb are the same 
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FIG. 10: Effective single-particle energies calculated for nuclei 
with an N=16 core using two Skyrme functionals that imple¬ 
ments tensor forces (Skxtb) or not (Skx). The down bending 
of the unbound IP 3/2 orbital observed when reaching the drip 
line is due to the halo like structure of the nuclei. 


TABLE IV: Effective single-particle energies for one neutron 
outside of closed shells for 24 O and 30 Si. Results are given for 
the configuration-interaction (Cl) and energy-density func¬ 
tions (EDF) models with the interactions discussed in the 
text. A is the shell gap at N=20. 


core 

Model 

e(04/2) 

(MeV) 

40 / 7 / 2 ) 

(MeV) 

4 IP 3 / 2 ) 

(MeV) 

A 

(MeV) 

0 

-7 

N 

EDF Skx 

-0.14 

5.14 

3.61 

3.75 


EDF Skxtb 

0.99 

4.13 

3.41 

2.42 


Cl WBA-M 

0.66 

4.88 

2.46 

1.54 

^sT 

EDF Skx 

-8.06 

-3.25 

-0.43 

4.81 


EDF Skxtb 

-8.06 

-3.25 

-0.43 

4.81 


Cl WBA-M 

-8.84 

-1.00 

-1.81 

7.03 


for Z=14 since the Skxtb tensor interaction for protons 
and neutrons cancel for the doubly jj closed shells in 28 Si. 
From Z=14 to Z =8 protons are removed from the 0 ( 4/2 
orbital. The tensor interaction between the 04/2 pro¬ 
tons and the 04/2 neutrons increases to its maximum 
for Z= 8 . This increases the ESPE of the 0 * 4/2 orbital 
and makes it unbound at Z= 8 . It also decreases the 
ESPE of the O/ 7/2 orbital. Thus, in the EDF model the 
N=20 shell gap at Z =8 is reduced due to a combination 
of the tensor-interaction effect for the 0(4 / 2 — O/ 7/2 ESPE 
spacing and the weak binding effect for the 0 ( 4/ 2 — IP 3/2 
ESPE spacing. 

The ESPE have been evaluated for the WBA-M inter¬ 
action at Z =8 and Z=14 by restricting the neutron core 
to have a (O 4 / 2 ) 6 ) (lsi/ 2) 2 neutron configuration. The 


Cl and EDF (Skxtb) results are compared in Table IV 
Results are relatively close for 24 0. But the Cl shell gap 
for 30 Si is larger than with EDF, explaining the earlier 


observation that states in 28 Na are calculated at too high 
energy (see Table IV). 


V. SUMMARY 


The spectroscopy of 28 Na has been investigated by 
means of the /3-decay of 28 Ne at GANIL/LISE and the 
in-beam 7 spectroscopy through the fragmentation of 
31,32 Mg beams at NSCL/S800. New positive parity states 
with J 7 r =3 + and 4 + are proposed at 691 and 968 keV, re¬ 
spectively, while new negative parity states are proposed 
at 1233, 1353, 1636 and 2493 keV with likely spin assign¬ 
ments J 7 r =2 — , J=3 _ , 4“ and 5“ respectively. Other neg¬ 
ative parity states are tentatively proposed at 1481, 1749 
and 1792 keV, the spin and parity assignments of which 
is more uncertain. Using these complementary methods 
all components belonging to the multiplet of states J’ r =l- 
4 + arising from the proton-neutron 0 d 5 / 2 - 0 d 3/ 2 coupling 
have been discovered. With the recent studies of the 
same multiplet of states in the 26 F and 30 A1 isotones, the 
evolution of the binding energy of the J 7 r =l-4 + multi¬ 
plet has been compared to shell model predictions using 
the USDA interaction. While the relative energies of the 
J 7 r =l-4 + states are well reproduced with the USDA in¬ 
teraction in the N=17 isotones, a systematic global shift 
in binding energy by about 500 keV is observed when 
moving from the valley of stability in 30 A1 to the drip 
line in 26 F. The origin of this change may arise from a 
change in the proton-neutron 0 d 5 / 2 - 0 d 3/ 2 effective inter¬ 
action when exploring large proton to neutron binding 
energy asymmetry. Other possible reasons are proposed 
as well in the text. 

The presence of a multiplet of negative parity states 
J 7 r =l-4 _ around 1.5 MeV, likely arising from the 7 ^ 5 / 2 - 
z/p 3 / 2 coupling, as well as a tentative observation of 
a J=5 _ state around 2.5 MeV, likely arising from the 
7 rd 5 / 2 -Kf 7/2 coupling, confirm the collapse of the N=20 
gap and the inversion between the neutron f 7/ 2 and P 3/2 
levels when removing protons in the d 5 / 2 orbital toward 
the drip line. These states are globally more bound 
than calculated by about 300 keV, a feature that may 
be due the high dimensionality of the basis or/and to an 
slightly wrong determination of the effective interactions. 
These features have been discussed in the framework of 
Shell Model and EDF calculations, leading to the con¬ 
clusions that no bound negative parity state would be 
present in 26 F and that the 3/2 + ground state and 3/2“ 
first excited states would be separated by only about 
500 keV in 25 O. It is important in the future to con¬ 
firm the spin assignments as well as the structure of the 
proposed negative parity states in 28 Na using for instance 
the 27 Na(d,p) 28 Na reaction. 
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